Controlling and Detecting Spin Correlations of Ultracold Atoms in Optical Lattices 
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We report on the controlled creation of a valence bond state of delocalized effective-spin singlet 
and triplet dimers by means of a Dichromatic optical superlattice. We demonstrate a coherent 
coupling between the singlet and triplet states and show how the superlattice can be employed to 
measure the singlet-fraction employing a spin blockade effect. Our method provides a reliable way 
to detect and control nearest-neighbor spin correlations in many-body systems of ultracold atoms. 
Being able to measure these correlations is an important ingredient to study quantum magnetism in 
optical lattices. We furthermore employ a SWAP operation between atoms being part of different 
triplets, thus effectively increasing their bond-length. Such SWAP operation provides an important 
step towards the massively parallel creation of a multi-particle entangled state in the lattice. 

PACS numbers: 03.67.Pp, 67.85.Jh, 75.10.Kt 



Strong correlations in quantum many-body systems 
are a cornerstone of modern condensed-matter physics. 
They underlie the Mott insulator (MI) state of elec- 
trons in the cuprates, which feature an extremely rich 
phase diagram p] and exhibit high-T c superconductiv- 
ity upon doping [3J. The experimental realization of 
the MI with ultracold bosonic [3] and - more recently - 
fcrmionic atoms 4J in optical lattices has demonstrated 
the prospect of these system to address fundamental 
condensed-matter problems [15! . One crucial require- 
ment for the study of quantum magnetism with ultra- 
cold atoms is a sensitive probe of spin correlations that 
characterize magnetic phases and can be employed to de- 
termine the entropy of the system 6]. Furthermore, the 
controlled manipulation of nearest-neighbor spin correla- 
tions might enable to create low-entropy spin-correlated 
states that can e.g. be adiabaticaly connected to coupled 
dimer states or an antiferromagnetically ordered state, 
circumventing cooling problems [7] . 

In this paper, we demonstrate both the control and 
the detection of short-range spin-correlations with ul- 
tracold bosons in an optical superlattice. We create a 
three-dimensional array of effective-spin triplet pairs and 
induce coherent singlet-triplet oscillations (STO) on the 
bonds. We make use of the different parity of the sin- 
glet and triplet wavefunctions to distinguish the two after 
merging pairs of sites [5] . The underlying mechanism is 
the singlet-triplet blockade known from double-electron 
quantum dots [9]. The detection procedure can be ap- 
plied to directly measure the amount of nearest-neighbor 
singlet and triplet correlations in a two-species MI of neu- 
tral atoms making it a valuable method to measure the 
entropy, as for example the singlet density in a fermionic 
MI at half filling is expected to increase with decreas- 
ing temperatures |10j . We furthermore employ a mas- 
sively parallel SWAP gate between neighboring triplet 
pairs which provides an important step towards the gen- 
eration of multi-particle entangled states with possible 



application in measurement based quantum computation 

Our experimental setup consists of a three dimen- 
sional optical superlattice filled with 87 Rb atoms |12) . 
Along the x-direction, the superlattice is formed by two 
collincar retro-reflected laser beams of wavelengths X xs = 
765 nm (short lattice) and X x i = 1530 nm (long lattice). 
The relative phase cj> of the two standing waves can be 
adjusted freely to yield a potential of the form V(x) = 
V xs cos(4k x x) — V x icos(2k x x + 4>) with k x = 2n/X x i 
and where the lattice depths V xs and V x i are control- 
lable by the intensity of the individual laser beams (see 
Fig. la). Transverse monochromatic lattices with wave- 
lengths \ y z = 843 nm complete the three-dimensional 
array of double- wells (DWs). Furthermore, a magnetic 
field gradient d x B x can be applied to create a potential 
bias in the DWs which depends on the Zeeman state of 
the atoms. 

We now consider the situation of symmetric DWs (<fi — 
0), occupied by two atoms in two internal states labeled 
|t) and If the largest energy scale in the problem 
is given by the trap frequencies of the individual wells, 
the system can be described by a two-site Bose-Hubbard 
type model 

k =- J E ( & lL^R + h - C ) 

<r=T4 

+U/2 (n L (h L - 1) + h R {n R - 1)) 
+A m /2 (n TL - n iL - hf R + n± R ) , (1) 

with the tunnel coupling J, the on-site interaction energy 
U (see Fig. la) and the state-dependent bias of magni- 
tude A m which reflects the magnetic field gradient. The 
operator a a L(R) annihilates a particle in the spin state a 

localized in the left (right) well, h a> i,{R) = ®lrL(R)®aL(R) 
counts the number of particles per state and well and 

For strong repulsive interactions (U 3> J), the ground- 
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FIG. 1. a Schematic drawing of the DWs formed by the super- 
lattice, b Eigenstates of Eq. (2) for two particles versus A m . 
At zero gradient, the states |f, t) and ||, f) couple to yield the 
triplet state [to) and the singlet |s). c Relative eigenenergies 
for two 87 Rb atoms in the \F = = ±1} Zeeman states 

in the DW potential versus barrier height V xs - The notation 
[vi , v%] for V XB — > refers to the vibrational quantum numbers 
for the first and second particle. The insets show the overlap 



state energy Eq as a function of V xs as obtained from 
an extended two-site Hubbard model [19]. As the bar- 
rier is ramped down, the spin-symmetry of the state 
is preserved. While the triplet states are adiabatically 
connected to the two-particle vibrational groundstate 
[0,0]|to) of the underlying long-lattice well, the singlet 
state requires one particle to occupy the first excited or- 
bital, thus [0, l]|s). This can be seen as an analog to 
the spin-blockade in electronic quantum dots [§]. To ful- 
fill the bosonic statistics, the spatial symmetry of the 
two-body wavefunction has to match the spin-symmetry 
|20) . In consequence, merging pairs of sites and subse- 
quently measuring the number of band excitations can 
be used to distinguish triplet and singlet correlations in 
a two-species MI. The on-site exchange splitting between 
[0,1] |s) and [0,1] \t ) could previously be measured via 
on-site exchange oscillations 21 . The ability to coher- 
ently couple \to) to |s) by means of the gradient term A m 
further provides the possibility to distinguish the states 
|t, |) and ||, f) from \t ) and |s). 

We begin our experiments by loading a BEC of about 
9 x 10 4 87 Rb-atoms in the \F = l,mp = —1) Zeeman 
state from a magnetic trap with high offset field into the 



[c t0 | of the groundstate with the triplet \t ) and the amount 3D ^ ^ y = ^ ™ s &nd y 35 E v 

of double-occupancy pdbi- 

state manifold is characterized by an occupation of one 
particle per site [13] . Virtual tunneling processes via 
higher-energy states mediate an effective superexchange 
coupling in this subspace [T51 HE]- The corresponding 
effective Hamiltonian can be derived from Eq. (1) [17] : 



H, 



off 



J ex 

(Sl • S R + S R ■ Si) + A m (S , £ - S R ) , (2) 



with the (ferromagnetic) superexchange coupling J ex — 
2J 2 /U > and the effective spin-1/2 operators Sf — 

(4*%i + 4i a ti)/ 2 ' S i = &\Ai ~ 4i a ti)/ 2i and ^ = 
(h^i — n±i)/2. This simple two-site model is well known 
within the framework of electrons in double-dots [S] . The 
groundstate of the Hamiltonian Eq. (2) for J ex > 
and A m = is the three-fold degenerate effective-spin 
triplet which consists of |t_i) = ||,|), |t + i) = |t>t) 
and |t > = (| t4) + I t4))/-\/2. The singlet state 
|s) = (|t, l> - |t4))/\/2 is higher in energy by 2J ex 
(see Fig. lb). For A m ^> J ex , the degeneracy of the 
states |t, |) and ||, t) is lifted and they become the eigen- 
states replacing |s) and \t ). A coherent coupling of \s) 
and \to) is realized by rapidly changing A m from zero to 
A m ^> J ex , projecting onto these new eigenstates. The 
subsequent phase-evolution describes singlet-triplet oscil- 
lations with a frequency vsto — 2A m //i. Previously, the 
reverse situation of a coherent superexchange coupling of 
the states | f, t) and | t, t) a t A m = has been realized 
in the same system [T5] . 

When the barrier is removed adiabatically (V xs — > 0), 
the singlet and triplet states in the DW are transferred 
into the two-particle eigenstates of the underlying long- 
lattice well. In Fig. lc, we plot the eigenenergies of 
the two particles in the DW with respect to the ground- 



|22j . The resulting state is a MI with at most one atom 
per site. By raising the long lattice to V x i — 40 E* 1 
and removing the short lattice, we merge pairs of lattice 
sites to yield constructed atom pairs in the vibrational 
groundstate in the long-lattice wells [19] . We use a radio- 
frequency rapid adiabatic passage (rfRAP) to transfer 
the atoms into the |1,0) Zeeman state and afterwards 
switch off the magnetic trap, maintaining a homogeneous 
offset field of ~ 1.2 G. Atom pairs in 1 1,0; 1,0) are sub- 
sequently transferred into |1, +1; 1, — 1) triplet pairs by 
means of spin-changing collisions (SCC) 23 . The two 
Zeeman states |1, +1) and |1, —1) serve as the two effec- 
tive spin states |f) and ||) [18]. 

We split the long-lattice wells into symmetric DWs by 
ramping up the short lattice to 40 E^ s again within 10 ms 
which creates an array of delocalized triplets \to). When 
the superexchange coupling is fully suppressed, we switch 
on a magnetic field gradient of variable strength for a 
holdtime <hoid- We detect the emerging STO by ramping 
down the barrier and employing a band-mapping tech- 
nique [23]. A pulsed magnetic field gradient at the be- 
ginning of the expansion separates the different Zeeman 
states. In order to avoid spurious signals due to imperfec- 
tions of the rfRAP and the SCC, the band-excitations are 
measured for the \F = l,mp — +1) Zeeman state alone 
- the only one not present before the SCC. The loading of 
the array of triplets and the subsequent conversion into 
singlets yields the analogue of a valence bond solid (VBS) 
state with bosons [U [5S] [55] . This state is characterized 
by a bond order along the superlattice direction with the 
period of the long lattice. 

In Fig. 2a, we plot the measured fraction of band exci- 
tations after merging n exc (ihoid) for the \F = l,mp = 1) 
Zeeman state for a gradient of about lOG/cm. We ob- 
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FIG. 2. a Plot of the relative population in the vibrationally 
excited band n eX c(ihoid) for a gradient current of 80 A (circles) 
together with the fit of a sine- wave (solid line). The insets 
show typical absorption images obtained by the detection se- 
quence, b Fitted oscillation frequencies vsto versus gradient 
current (filled circles) together with an independent measure- 
ment of the single-particle shift v$p = A m //i (open circles). 
The solid line is a linear fit to vsto which translates into a 
slope of 0.122(1) G/(cm A) and an offset of -0.12(3) G/cm. 

serve an amplitude of the STO of about 30%. The re- 
duction from the ideal value of 50% stems from a resid- 
ual magnetic field gradient present during the splitting, 
the finite lifetime of the triplet and singlet state and im- 
perfections of the detection method. About 4% of band- 
excitations created during the loading procedure give rise 
to an additional offset. The phase-shift of the oscilla- 
tions is a result of the finite switch-on time of the gra- 
dient. Fig. 2b shows the measured oscillation frequency 
vsto which depends linearly on the current applied to 
the coil producing the gradient. We independently mea- 
sure the energy bias hv$p between the left and the right- 
well for a single particle versus gradient strength. This 
is done by coherently splitting long-lattice sites with sin- 
gle atoms in |1) and recording the emerging interference 
patterns after a short holdtime ~ 1 ms and time-of-flight 
[IB] . We find the STO frequency to be 2v S p (see Fig. 2b) 
which confirms the two-body nature of the measured ef- 
fect. From a linear fit of ^sto we find a gradient offset 
of ~ 120mG/cm which might stem from inhomogencitics 
of the 1.2 G offset field and static magnetic field sources 
close to the experiments. An alternative explanation is a 
relative angle between the polarizations of the two lattice 
beams of the short lattice of < 1 mrad leading to a weak 
spin-dependency of the superlattice [27] . 

For all gradient strengths, we find a damping time of 
the oscillations of about r = 40 — 50 ms. This damping is 
mainly caused by the decay of the triplet and the singlet 
state. We measure the lifetime of the triplet by merging 
the DWs after the holdtime ihoid without applying the 
magnetic field gradient. The same is done for the singlet 
state, where only a short gradient pulse after the split- 
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FIG. 3. Measured lifetime of the triplet (dark blue) and the 
singlet state (light blue) versus holdtime in the split DWs. 
The blue solid lines are the result of a simultaneous exponen- 
tial fit yielding a lifetime of 43(1) ms. The black data points 
show STO taken with the same sequence and an additional 
gradient field of about 10 G/cm during the holdtime. 



ting is employed to convert |trj) into |s). In Fig. 3, we 
plot n exc for the two measurements as a function of ihoid 
together with a STO trace taken at a gradient current 
of 80 A. The lifetime measurement provides an envelope 
to the oscillation data. A simultaneous fit of the triplet- 
and singlet-lifetime yields r — 43(1) ms which matches 
the damping-time of the STO. This lifetime gradually de- 
creases when increasing the short-lattice and transverse- 
lattice depths. Neither spontaneous scattering of lattice 
photons (r^ 1 > 500 ms), nor tunneling of individual par- 
ticles can explain this lifetime. It is most likely limited by 
weak spin dependencies of the external confinement due 
to imperfections in the polarization of the lattice beams 
[T§] . A larger detuning of the lattice beams would signif- 
icantly reduce the sensitivity on these imperfections. 

In addition to the formation and detection of the delo- 
calized triplets, the superlattice also offers unique possi- 
bilities for further manipulation (Fig. 4a), i.e. aimed at 
generating multi-particle entangled states or to increase 
the spatial extension of the entangled spin pairs [28] . Af- 
ter having created the array of | to) bonds on neighbor- 
ing sites, we remove the long-lattice, shift its phase to 
(j> = 7r and ramp it high again, thus combining pairs of 
sites which belong to different triplet-bonds to a DW. 
By lowering the short lattice to 12E* S within 200 /is, we 
switch on a superexchange coupling between these sites 
with J ex = J/3 = h x 360 Hz [TS]. We ramp the bar- 
rier high again after half a superexchange period which 
realizes a SWAP operation for the coupled spins. As a 
result, the triplets are now delocalized over three lattice 
spacings rather than a single one before [35]. We again 
use the magnetic field gradient to induce the STO and 
subsequently repeat the SWAP operation, in order to re- 
store the original bond length of the entangled spin pairs. 
Finally, the same combination of merging and band map- 
ping is used as before to reveal the STO. The measured 
fraction n ea;c (ihoid) is plotted in Fig. 4b together with a 
trace recorded with the same sequence but without low- 
ering the short lattice to induce the two SWAPs. We 
find an oscillation frequency three-times higher when the 
SWAP operations are carried out. This is explained by 
the linear dependence of A m on the distance of the par- 
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FIG. 4. a Scheme of the SWAP-sequence: The spin pairs 
(1) are split into triplets (2). A superexchange for half a 
period SWAPs the spins of neighboring triplets (3) before a 
magnetic field gradient is used to drive STO (4). A second 
SWAP (5) brings the spins back to their original position 
(6). b Measured excited fraction after the final merging with 
(blue circles) and without (black circles) applying the SWAP 
gates. The oscillations with the SWAP are faster by a factor 
of 2.91(10), confirming the successful stretching of the spin- 
pairs. 

tides for a given gradient strength. The corresponding 
STO amplitude is reduced by about 40% which in parts 
stems from unsuccessful SWAPs where there is a hole ad- 
jacent to a triplet bond. These holes are detected as ad- 
ditional excitations, yielding a higher offset of the STO. 

The successful demonstration of the SWAP operation 
constitutes an important step to further entangle neigh- 
boring triplet pairs which can be achieved by a VSWAP 
operation. Since this manipulation is carried out between 
all neighboring singlets simultaneously, the result after a 
single step would be an entangled chain of neutral atoms. 
With a second superlattice along a perpendicular direc- 
tion, this entanglement can be extended on a 2D plane 
in the same manor. The resulting state is maximally en- 
tangled in 2D and therefore would be valuable for mea- 
surement based quantum computation. 

In conclusion, we have demonstrated the controlled 
loading of an array of nearest-neighbor S z = triplets 
and the conversion into singlets by means of an optical 
superlattice and a magnetic field gradient. The array of 
singlet or triplet dimers can be seen as a valence bond 
solid with the bonds alternating along the x-axis. We 
note that the same method can be applied to create a cou- 
pled dimer antiferromagnet of fermions |26j from an ini- 
tial low-entropy band-insulating state. We have also es- 
tablished a sensitive method to measure nearest-neighbor 
spin correlations to reveal gradient driven singlet-triplet 
oscillations. Finally we have demonstrated the controlled 
stretching of the triplet pairs by a massively parallel 
superexchange SWAP gate between neighboring bonds. 
Bringing distant spins together in a DW by SWAP opera- 
tions is a way to measure longer-range spin-correlations. 
The implementation of a %/SWAP gate instead can be 
used to entangle chains - and further 2D arrays - of neu- 
tral atoms. 
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APPENDIX I - EXTENDED BOSE-HUBBARD 
MODEL 

We calculate the eigenenergies in the double-wells 
as a function of V xs (Fig. la) using an extension of the 
Bose-Hubbard model Eq. (1). Since the overlap of the 
left- and right-localized wavefunctions becomes signifi- 
cant for low barriers - thus small values of V xs - a direct 
nearest-neighbor interaction needs to be taken into ac- 
count, as well as a density-dependent modification of the 
tunnel coupling |18j . The corresponding extended Bose- 
Hubbard Hamiltonian reads 

Hebhm = H - A J ^2 (haL + n aR ) {a\, L a a , R + h.c.) 

1 .+ . . 

+ 2 a aL a a'L a a'Ra aR 
+ 2 a <jR a a>R a <?'L a aL ) , 

with Ulr. = g x J w 2 l {t)w r (r) dr and AJ = —g x 
/ w^(r)w R (r) dr. The functions wl,r(tc) are the left- 
and right-localized groundstate wavefunctions functions 
in the double-well. The prefactor to the integrals is 
g = 4ir 2 h 2 as /wRb with ag ~ 5.3 nm being the s-wave 
scattering length and raRb the mass of the 87 Rb atoms. 
All parameters of the extendend Bose-Hubbard model 
Eq. (3) are obtained from a single-particle bandstructure 
calculation. 



APPENDIX II - LOADING OF ATOM PAIRS IN 

\F = l,m F = 0) 

During the rampup of the lattice 3D optical lattice to 
load the MI, we apply a shallow long-lattice potential of 
V x i — 1 £y with a phase 4> — 0.2. Therefore, possible 
defects are loaded mostly to "left" sides of the double- 
wells. When merging the pairs of atoms, we apply a small 
bias A < U {V xl = 4QE?, 4> ~ 0.05), to avoid band- 
excitations from these singly occupied double-wells. As 
a result of the loading sequence, we obtain ~ 60% of the 
atoms loaded into pairs. 

The offset-field of the magnetic trap from which the 
optical lattice is loaded is about 153 G. After loading 
the MI, we transfer all atoms from \F — 1,to^ = — 1) 
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to 1 1,0) by coupling them with a radio-frequency sig- 
nal of 108.6 MHz and sweeping the magnetic field across 
the resonance at ~ 152.5 G. We optimize the mag- 
netic field ramp to obtain no population in |1,1) and 
achieve a transfer efficiency of more than 96%. Due to 
the quadratic Zeeman shift, the transition |1,0) -H> |1, 1) 
is detuned by about 3 MHz. After the rfRAP, the high- 
offset trap is switched off, leaving a homogeneous offset 
field of - 1.2 G. 



APPENDIX III - LIFETIME OF THE SINGLETS 
AND TRIPLETS 



where, k VtZ = 2n /Xy, z is the transverse-lattice wavenum- 
ber. If the polarization of the beams has a non-zero circu- 
lar component with respect to the quantization axis (i.e. 
the magnetic field), the transverse-lattice depth V VtZ will 
depend on the internal state of the atoms [57] ■ There- 
fore, the harmonic confinement along the x-direction will 
be spin-dependent which translates into a spatially inho- 
mogeneous distribution of A m . Given our experimental 
parameters and the measured decrease of r with V y . z , 
we estimate a relative difference of the transverse-lattice 
depths for the states |f) and \\) of ~ 0.6%. 



In order to explain the lifetime r of the singlets and 
triplets observed the experiment (see Fig. 3) we have 
measured the dependence of r on the transverse-lattice 
and short-lattice depths. In both cases, we find a grad- 
ual decrease of r with increasing lattice depth. While 
the decrease in lifetime with short-lattice depth is com- 
patible with the increase in the spontaneous scatter- 
ing rate of lattice-photons, the dependence on V ViZ is 
much stronger. We explain this behavior by assuming 
a small spin-dependency of the transverse lattices caus- 
ing a spin-dependent confinement in the x-direction due 
to the Gaussian shape of the lattice beams. 

The red-detuned transverse lattice beams have a waist 
of Wy^ z ~ 150 /im. The confinement due to their Gaussian 
profile can be approximated by a harmonic potential with 
trap frequency 



The quantization axis in our experiments is mainly 
oriented along the x-direction, thus minimizing spin- 
dependencies in the transverse lattices. Small angles a 
between the quantization axis and the ^-direction, how- 
ever, are possible. Furthermore, birefringence in the 
optics passed by the initially linearly polarized lattice 
beams can yield the required ellipticity of the polariza- 
tion. For a ~ 10, an ellipticity of e ~ 0.1 in the trans- 
verse lattice beams leads to the aforementioned spin- 
dependency of the lattice depth. These values are slightly 
reduced when also considering the dependency of the 
spontaneous scattering rate in the lattice on V y ^ z . We 
finally note, that the dephasing due to a spin-dependent 
external confinement is stronger for longer-ranged spin- 
correlations. 
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